Intraband absorption in finite, inhomogeneous quantum dot stacks for 
intermediate band solar cells: limitations and optimization 
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We present a theoretical analysis of intraband optical transitions from the intermediate pseudo-band of 
confined states to the conduction band in a finite, inhomogeneous stack of self-assembled semiconductor 
quantum dots. We show that the absorption spectrum already for a few coupled dots differs from that of a 
single dot, which leads to an enhancement of the overall transition rate under solar illumination. The decisive 
role of the strength of inter-dot coupling for the stability of this enhancement against QD stack inhomogencity 
and temperature is revealed. 
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One of the ways to improve the efficiency of solar cells 
is to introduce an intermediate band in the energy spec- 
trum of a photovoltaic structure 1,2 . In this way, electrons 
can be sequentially promoted from the valence band to 
the intermediate band and then to the conduction band 
by absorbing photons with energies below the band gap 
which arc not converted into useful electrochemical en- 
ergy in a standard structure. As an implementation of 
this concept, a stack of quantum dots (QDs) in the intrin- 
sic region of a p-i-n junction solar cell has been proposed 3 . 
This idea has indeed gained some experimental support 
in recent years 4-12 . Quantum-dot-embedded p-i-n so- 
lar cells show higher quantum efficiency in near infrared 
range but their overall efficiency still is lover than the 
efficiency of similar devices without QDs 4-12 . 

On the theory side, models involving a single QD were 
formulated to describe the kinetics of transitions from 
and into the intermediate levels 13,14 . On the other hand, 
modeling of the electron states and optical absorption in 
chains and arrays of QDs has been mostly limited to in- 
finite, periodic superlattices of identical dots 15-20 . As we 
have shown recently 21 , enhanced absorption can appear 
also in finite chains of non-identical QDs but it is sup- 
pressed if the inhomogencity of the QD chain (leading 
to non-identical electron ground state energies in the in- 
dividual QDs) becomes too large. Since the actual QD 
chains arc always finite (usually built of several to a few 
tens of QDs) 4-12 and unavoidably inhomogeneous it is of 
large practical importance for the optimal design of in- 
termediate band photovoltaic devices to extend the the- 
oretical analysis to such more realistic structures. 

In this paper, we study the intraband optical absorp- 
tion associated with the electron transition from the 
states confined in a finite stack of quantum dots to the 
conduction band. We propose a relatively simple and 
computation-effective model which, however, includes all 
the essential features of the system, in particular the 
inhomogencity of the energetic parameters of the dots 
forming the stack and the coupling between them. We 
show that the enhanced absorption features appear al- 
ready for a few dots and lead to enhanced transition 



rate from the pseudo-band of confined states to the bulk 
continuum. While this effect is relatively stable with re- 
spect to inhomogeneity, it turns out to be associated with 
the transitions from the ground state of the stack and is 
washed out as soon as the temperature becomes compa- 
rable to the width of the pseudo-band of confined states. 
Our modeling results suggest that this detrimental tem- 
perature effect can be to a large extent overcome by in- 
creasing the tunnel coupling between the dots. 

We are interested in the transitions from a single- 
electron bound state \v) (with a wave function ^ y (r) 
and energy E u ) to a continuum state |fc) (with a wave 
function ^ k (r) and energy E k = h 2 k 2 /(2m*)). We use 
the dipole Hamiltonian for the coupling to incident light, 
Hint = e l {£o£oo)r ■ ££o cos(q • r — uit), where £n is the vac- 
uum permittivity, is the high-frequency dielectric con- 
stant of the semiconductor, e is the elementary charge, 
£0 is the amplitude of the electric field of the electro- 
magnetic wave, £ is a unit vector defining its polariza- 
tion, q is its wave vector (inside the dielectric medium), 
and to = cq/n T is its frequency, where n r is the refrac- 
tive index of the semiconductor. Using the Fermi golden 
rule, we find the transition rate for the illumination with 
monochromatic radiation 
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a„ k (q)\ 2 S(hu + E„-E k ), (1) 



where U = £o£oo-Eo/2 is the energy density of the elec- 
tromagnetic wave and the transition matrix elements are 
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Under illumination by a broad band radiation, the cor- 
responding transition rate per unit frequency interval is 
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where u(lj) is spectral density of radiation. 

In our model, we assume that the confined electron 
states \v) are superpositions of the ground states \n) 
confined in the individual dots (where n numbers the 
dots). For simplicity, we assume that each of these sin- 
gle dot states has an identical wave function, ip n [ r ) = 
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ipo(r — Rn), where i?„ is the position of the nth dot (we 
assume that the dots are stacked along the growth direc- 
tion z). The ground state electron energies in the dots, 
e„, may differ. The states \n) are coupled by nearest 
neighbor couplings. The eigenstates \v) and the corre- 
sponding energies E v are thus obtained as the eigenstates 
of the effective chain Hamiltonian (assuming a single con- 
fined state in each dot) 21 , 



H = 



N-l 

E 

n=0 



JV-2 



e„|n)H + t J2 d"X" + M+ h.c), (4) 
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where t is the coupling constant and we assume the over- 
lap between the wave functions localized in different dots 
to be negligible, so that (n\n') = 5 nn >. The inhomogene- 
ity of the QD stack is taken into account by choosing the 
energies e n from the Gaussian distribution with the mean 
E and variance a 2 . 

We assume that the wave function for the elec- 
tron in the nth dot has the Gaussian form, ty n (r) = 
AAexp{-(l/2)[(x 2 + y 2 )/l 2 + (z - z„) 2 // 2 ]}, where z n = 
nD is the position of the nth dot, I, l z are the extensions 
of the wave function in the xy plane and along z, respec- 
tively, and AT is the normalization constant. For the bulk 
electron states, we assume plane waves 22 orthogonalized 
to the localized states, as previously proposed for calcu- 
lating carrier capture rates 23-25 . These states are labeled 
by the wave vector k describing the plane wave far away 
from the QD structure. Thus, we have 
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where Nk is the appropriate normalization constant, we 
assume normalization in a box of volume V with peri- 
odic boundary conditions, and the orthogonalization co- 
efficients 7fc„ are given by 
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In all our simulations we focus on an InAs/GaAs struc- 
ture (not optimal for a solar cell 1 but important for the 
current laboratory scale investigations 2 ). Wc assume the 
wave function confinement sizes I — 4.5 nm, l z = 1.0 
nm (see Ref. 21 for the discussion of the dependence on 
these parameters), E = —250 meV. Based on earlier k-p 
calculations 26 , the tunnel coupling t is given by the for- 
mula ln(t/to) = —/CD, where JC = 0.59 nm" , to = 0.79 
eV (t = —3.9 and —41.3 meV for D = 9 and 5 nm, re- 
spectively) . 

In Fig. 1(a), we show the intraband absorption spectra 
of short chains of identical QDs compared to the spec- 
trum of a single dot for light incident along the chain (in 
the z direction) and perpendicular to the chain. Both 
spectra are calculated assuming that only the lowest con- 
fined state is considerably occupied (which corresponds 
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FIG. 1. (a) Low temperature intraband absorption spectra 
for a chain of identical dots illuminated by unpolarized light 
incident along and perpendicular to the stacking direction z 
(as indicated by e z and e x , respectively), (b) The dependence 
of the spectrum on the angle of incidence for 10 QDs. 
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FIG. 2. Low temperature intraband absorption spectrum for 
a chain of non-identical dots illuminated by unpolarized light 
incident along the stacking direction z. (a, b) As a function of 
the degree of inhomogeneity; (c, d) as a function of the chain 
length. Each spectrum is an average of 1000 realization with 
randomly chosen values of e n . 



to low temperatures) and for unpolarized light (for per- 
pendicular incidence, the contribution from light polar- 
ized along z is small). In both cases, the absorption spec- 
trum for a single dot shows only a single maximum with 
a long tail on the high energy side but already for a few 
dots it develops a series of additional maxima. As can 
be seen in Fig. 1(b), the spectra for a chain of 10 QDs 
smoothly evolve between the two limiting geometries as 
the incidence angle is tilted from normal to in-plane in- 
cidence. 

The additional absorption peaks that can be seen in 
Fig. 1 result from interference effects which are due to 
derealization of the electron state and lead to preferred 
transitions to states with k z = 2irj/D, where j is an 
integer 21 . Disorder, which in our case has the form of 
inhomogeneity of the "on-site" energies e„, destroys the 
coherently delocalized electron states and leads to local- 
ization. As a result, the additional peaks are suppressed 
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FIG. 3. Intraband absorption spectrum for a chain of 10 non- 
identical dots with a — 7.8 meV, illuminated by unpolarized 
light incident along the stacking direction z as a function of 
temperature for two values of the inter-dot coupling: (a, b) 
t — —3.9 meV (corresponding to D = 9 nm), (c, d) t = —41.3 
meV (corresponding to D — 5 nm). 



and the absorption spectrum becomes more similar to 
that of a single dot (see Fig. 2(a, b)). As we show in 
Fig. 2(c, d), while the spectra for QD stacks differ con- 
siderably from those for a single dot, the extension of the 
stack above 10 dots has little effect. 

For a realistic device, absorption at elevated temper- 
atures is relevant. At non-zero temperatures, the occu- 
pation of excited states in the QD-related pseudo-band 
is non-negligible, which leads to a reconstruction of the 
spectrum, as shown in Fig. 3. For chains with \t\ < kT, 
the structure of the spectrum is washed away with in- 
creasing temperature (Fig. 3(a, b)). However, as can 
be seen from Fig. 3(c, d) strong coupling broadens the 
pseudo-band beyond the thermal energies and stabilizes 
the chain absorption against temperature. 

In order to study the properties of the device under re- 
alistic working conditions we calculate also the electron 
intraband transition rate /?o (from the ground state) for a 
chain of non-identical QDs illuminated by black body ra- 
diation. The energy density of radiation with T = 5777 K 
was normalized to the incident flux of 1 kW/m 2 (one 
sun). In Fig. 4(a,b) we show the spectral distribution of 
the transition rates as a function of the energy of the ab- 
sorbed photon. Comparison with the absorption spectra 
shown in Fig. 1 shows that the qualitative form of these 
two spectral characteristics is similar but quantitatively 
the transition rate at higher photon energies is enhanced 
relative to that at lower energies. This is due to the very 
quickly growing spectral densities of the thermal photon 
field at higher energies (according to the the Planck law) 
which enhances the role of higher energy photons (and, 
in consequence, the overall absorption of energy in the 
corresponding spectral range). One should notice that 
inhomogencity (higher values of a) leads to a less struc- 
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FIG. 4. (a, b) Electron intraband transition rates under il- 
lumination by one sun thermal radiation as a function of the 
degree of inhomogeneity (a) and as a function of the num- 
ber of QDs in the chain (b). Each spectrum is an average 
of 1000 realization with randomly chosen values of e„. (c, d) 
Integrated single-electron transition rate under illumination 
by one sun thermal radiation as a function of the degree of 
inhomogeneity (c) and as a function of the number of QDs in 
the chain(d). 



tured spectrum with non-vanishing absorption over all 
the spectral range. 

In fig Fig. 4(c), we show the total transition rate, in- 
tegrated over photon energies. The result at 300 K in- 
cludes the thermal distribution of the initial states of 
the electron. As one can see, the total rate increases as 
the inhomogencity grows. This increase is due to the 
increasing probability that the electron will localize in 
a dot with much lower "on-site" energy as the QD en- 
ergy distribution broadens. This enhances the transition 
rate because of the growing spectral density of radiation 
for higher energy differences between the initial and fi- 
nal states. As a function of the chain length (Fig. 4(d)), 
the total transition rate increases slightly (up to several 
%) for strongly coupled chains at low temperatures and 
then saturates for N > 10. This increase is mostly due 
to high-frequency contributions corresponding to the en- 
hanced high-energy absorption features in a QD stack (cf. 
Fig. 1). One can notice in Fig. 4(d) that the gain in the 
total transition rate for strongly coupled dots not only is 
much larger than for weaker coupling but it is also much 
less sensitive to temperature. 

In summary, we have proposed a model which describes 
electron states and intraband absorption spectra which 
correspond to the transitions from the QD pseudo-band 
to the conduction band. We have studied the absorp- 
tion spectra as a function of the degree of inhomogene- 
ity, the number of QDs in the chain, the inter-dot sepa- 
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ration, the temperature and the illumination conditions. 
We have shown that the absorption spectrum even of 
short QD chains (a few QDs) is dominated by interfer- 
ence effects that leads to the appearance of additional 
absorption maxima and does not considerably evolve fur- 
ther if the chain length exceeds approximately 10 dots. 
The overall transition rate under illumination with black 
body radiation increases by up to several per cent with 
the number of QDs but also saturates already for a few 
QDs in the chain. The structured absorption spectrum 
persists when the inhomogencity of the QD transition 
energies is increased up to the values comparable with 
the magnitude of the inter-dot coupling. Strong cou- 
pling between the dots is essential for maintaining the 
chain absorption features up to high temperatures. On 
practical side, our results show that already a stack of a 
few QDs manifests the absorption features characteristic 
of QD chains and that these features are stable against 
temperatures and energy inhomogeneities if the dots are 
sufficiently strongly coupled. However, a QD chain shows 
only a slightly increases absorption as compared to a sin- 
gle dot and the increase of absorption is dominated by 
high-energy photons, which may lead to competition with 
the interband absorption. Moreover, a trade-off has to be 
found between the absorption enhancement in strongly 
coupled chains and the need to reduce tunneling in order 
to avoid carrier escape from the intermediate bands 27 . 
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